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Abstract 

We discuss the relationship between the phase diagram of the Q = state Potts model, 
the arboreal gas model, and the supersphere sigma model S ' 2 = OSP(l/2)/OSP(0/2). We 
identify the Potts antiferromagnetic critical point with the critical point of the arboreal gas 
(at negative tree fugacity), and with a critical point of the sigma model. We show that the 
corresponding conformal theory on the square lattice has a non-linearly realized OSP(2/2) = 
SL{\/2) symmetry, and involves non-compact degrees of freedom, with a continuous spectrum 
of critical exponents. The role of global topological properties in the sigma model transition 
is discussed in terms of a generalized arboreal gas model. 

1 Introduction 

The question of possible fixed points and RG flows in 1 + 1 theories with supergroup symmetries 
plays a major role in the supersymmetric approach to phase transitions in non-interacting dis- 
ordered systems (such as the transition between plateaux in the integer Quantum Hall or spin 
Quantum Hall effects). Although some similarity with what happens in systems with ordinary 
group symmetries is expected, it has also become clear over the last few years that more com- 
plicated scenarios are at work. Among the features that play a key role in these scenarios are 
the facts that supergroup symmetries are not protected by the Mermin- Wagner theorem and can 
be spontaneously broken in two dimensions; that the lack of unitarity allows non Kac-Moody 



1 



Figure 1: The configurations of the arboreal gas model are forests of trees on the given graph; 
the union of trees jointly spans the lattice and each component tree carries a Boltzmann weight 
a. As an example, we show a forest on a portion of the square lattice. There are six component 
trees, including an isolated vertex. 

(logarithmic) conformal field theories with supergroup symmetry; etc. Recent references in this 
area are [H II El HE]- 

A systematic program to investigate these questions would be to tackle the continuum limit 
of spin chains with supergroup symmetries, either with compact or non-compact representations, 
with or without integrability. This turns out to be quite hard (for some recent work in this 
direction, see Refs. as profoundly new features appear, on top of, and maybe related to, 

technical difficulties, already in the first steps of this program. 

We discuss in this paper one of the simplest possible examples of a supergroup theory one may 
imagine, the supersphere sigma model OSP(l/2)/OSP(0/2). This model is deeply related to the 
arboreal gas model (see Fig.^), and the Q — > limit of the Potts model on the square lattice, and 
through these relations, can be studied in considerable detail. In particular, the relation with the 
Potts model provides a phase diagram for the sigma model, a confirmation that the perturbative 
results are identical with those of the 0(n),n = — 1 sigma model, and the identification of 'a' 
critical point of the sigma model with the antiferromagnetic transition in the Potts model, as 
well as the transition in the arboreal gas model. Surprisingly, the associated critical theory turns 
out to be of a type seldom if ever encountered before in lattice models: with c = —1, it involves 
a non-compact boson as well as a symplectic fermion, and is related to a theory with a larger, 
OSP{2/2), symmetry. 

The paper is organized as follows. In section 2, we review the relationship between the 
Q = Potts model, the arboreal gas model, and the sigma model. A short version of the results 
presented there has appeared in Ref. 0. Related results have also appeared in Ref. 0. Sections 
3 and 4 discuss the continuum limit of the model at the antiferromagnetic critical point using 
Bethe ansatz equations as well as techniques from conformal field theory. In section 5 we prove 
the remarkable identity between the Bethe equations for the Q = antiferro Potts model and the 
integrable osp{2/2) spin chain in the fundamental representation. Using general results about the 
phase diagrams of orthosymplectic chains provides independent confirmation of our continuum 
limit identification. Some conclusions are gathered in section 6 while the appendix discusses a 
few technical aspects of Bethe ansatz for the sZ(2/l) = osp(2/2) algebra. 

This paper is not only a step up the program of understanding statistical mechanics models 
with supergroup symmetries: first, the results will prove extremely useful in the analysis of the 
antiferromagnetic Potts model for Q generic, which we will discuss elsewhere [Sj; second, the 
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arboreal gas is a model of great interest from the point of view of transitions in gases of 'hard' 
objects at negative fugacity; finally, the relationship between the arboreal gas and the sigma 
model provides a way to investigate fundamental questions — like the role of global topological 
properties in sigma models in their broken symmetry phase — in a particularly simple context. 

2 The Q = Potts model and the sigma model 

The Potts model at Q = can only be well defined through the cluster expansion and a proper 
limiting procedure jlUj . This requires generally to approach coupling K = along lines e K — 1 = 
2, in which case the partition function as Q — ► can be rewritten as the partition function of 
an arboreal gas |H] 

oo 

^arboreal = E flPAr ( F p)' W 
p=l 

where N(F p ) is the number of ways of covering the lattice with p non-overlapping trees (a forest), 
such that each site belong to one and only one tree. 

The proof of this statement is straightforward. Write the partition function as 

^pott s = e n * KS ^ 

<Ti=l,...,<2 <ij> 

= £(e*-i)V' 

G 

« (^) S J2 aCl ' Cy Q Cy ( 2 ) 

G 

where G are subgraphs of the lattice made of B bonds and CI clusters (connected components), 
S is the total number of sites, Cy the number of cycles on the clusters. We have used Euler's 
relation (valid in the plane or on an annulus, a geometry we will consider below), S = B+Cl — Cy. 
In the last equation Q 0, e K — 1 ^. As Q — > 0, only graphs without cycles survive, giving 
the result Q after a renormalization of Z. This is true irrespective of the lattice. 

On the square lattice, the usual critical Potts model corresponds to approaching Q = along 
the selfdual curve e K — 1 = yfQ [TTj, i.e., by taking a — > 0, in which case one keeps a single 
spanning tree. The geometrical properties are then described in terms of a symplectic fermion 
theory [12111311111. 

Similarly, on the triangular lattice the selfdual curve reads (e K - l) 3 + 3(e^ - l) 2 = Q PH, 
which, for Q — > and K > 0, becomes e K — 1 ~ whence a — > and the equivalence with a 

spanning tree follows once again. These two lattices thus strongly suggest that the usual critical 
behaviour is universal, and governed by symplectic fermions. 

The case a > corresponds to a massive region in the Potts model phase diagram, and 
therefore the arboreal gas is massive. The case a < negative is more interesting. While the 
model is still massive for \a\ big enough, on the square lattice and for a* < a < 0, the model is 
massless, and flows to the symplectic fermion theory again. The special value a* = —4 is critical, 
and corresponds to the limit Q — > of the critical antiferromagnetic Potts model on the square 
lattice [HI [□]. 

On the triangular lattice there is strong numerical evidence for a massless regime when a* < 
o < 0, again characterized by a flow to the symplectic fermion theory 0. However, since 
a* = —5.705 ± 0.006 is only known numerically, it appears to be difficult to detect whether there 
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is distinct critical behaviour exactly at a*, and what is the exact nature of the transition to the 
massive regime a < a* [Hj- 

It turns out that, on any planar graph, the arboreal gas can be expressed in terms of an 
interacting symplectic fermion theory j7]. This comes from the remarkable identity 



/oo 
d7 ? id7 ?2 e a(Ql - Q2 W /,?2 = Y / a p N{F 1 
P =i 



(3) 



Here, each lattice site i carries a pair of fermions 771 (i) and rj2(i), dr]idr]2 is shorthand for the 
measure of integration over every site, and M is the discrete Laplacian with matrix elements 
rriij. Thus, for i ^ j, rriij equals minus the number of edges connecting i and j, and ma equals 
the number of neighbours of site i (excluding i itself, if the graph contains a loop). We have 
introduced the two objects 

Qi 
Q2 



(4) 
(5) 



(ij) 

where we recall S is the total number of sites on the lattice, and the sum in Q2 is taken over pairs 
of neighbouring sites. This fermionic reformulation is interesting for several reasons, in particular 
because it allows one to derive quick conclusions on the phase diagram of the arboreal gas. One 
observes first that the action in Eq. exhibits a non-linearly realized OSP(l/2) symmetry. 
Indeed, let us introduce on each vertex an auxiliary field (ft(i) subject to the constraint that 
4> 2 {i) + 207/1(2)7/2(2) = 1. Solving for <fi gives two solutions, of which we chose 0=1 — 0771772, 
that is the 'body' (part made of c- numbers) of <f> positive (which we denote by b(<f>) > 0). Using 
the basic rule that 5(4> 2 + 2ar\\r\2 — 1) = j$5{(f> + 0771772 — 1), together with ^ = e aVim , one can 
entirely eliminate the aQi term in the action and write the functional integral as 

•^arboreal 



2 N 



b(4>)>0 



x exp 



J|d0(i)d77i(7)d77 : 



1» 



(t) + 20771(^)772(7:) - 1] 



(ij) 



1,3 



(6) 



We now introduce the vector u = (0, 771,772) with the scalar product u(i) ■ u(j) = <j)(i)(j)(j) + 
0771 (7)772 (j) + 0771(^)772(7). By explicitly writing down the discrete Laplacian terms and regrouping 
a bit, one finds that the action can be rewritten as 



1 



S = J2(u(i)-u(j)-l) 



(7) 



(ij) 



The additional four-fermions term have disappeared from the action, which exhibits global 
OSP(l/2) symmetry (the group of transformations leaving the scalar product invariant). The 
generating function in fact looks exactly like the one of a discretized sigma model on the target 
manifold OSP(l/2)/OSP(0/2) = UOSP(l/2)/SU (2)1 

Although the parameter a appears in our definition of the scalar product, it can always be 
absorbed in a rescaling of the fermions, and thus 



O -^arboreal 



b(<j})>0 



[d0d77id7/2] 5((f) 2 + 2771772 - 1) exp 



1 



-5>(iMC;)-i) 



(ij) 



(8) 
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where now u(i) ■ u(j) = <j>(i)(j)(j) + f]i{i)f]2{j) + f]\ (j)V2 (*) ■ An infinitesimal OSP(l/2) transfor- 
mation reads 



5(f) = -SZirj! + 5&V2 
Srji = -5Z,2<t> + fi a Vi + ^ cr l2 

Sr] 2 = —5£i<j) + Sbrji - Sarj 2 (9) 

where 5£i,<5£2 are 'small' fermionic deformation parameters, and 5a, 5b, 5c are small bosonic 
parameters. In terms of the fermion variables, the symmetry is thus realized non-linearly: 

Srji = -<5£ 2 (1 - V1V2) + 8ar]i + 5crj 2 

5rj 2 = — — 771772) + Sbrji — Sarfe (10) 

We now take the continuum limit by introducing slowly varying fermionic fields r]i(x,y) and 
r]2 (x,y). The euclidean action reads after some simple manipulations (Boltzmann weight e~ s ) 

s = - J d 2 x [(d^f + 20^0^2] . (11) 

with coupling g = —2a. Note that no lattice spacing (say, do) appears in this action as the 
Oq from the expansion of the fields is compensated by the l/oo from the transformation of the 
discrete sum over lattice sites into an integral. 

We now go back to the issue of the sign of the field <f>. There are of course two solutions to the 
equation 4> 2 + 2r\\r\2 = 1. While we have focused on eft = 1 — 771772, the other choice <fi = —1 + 771772 
is also possible. It is therefore more correct to say we have mapped the Q — > Potts model 
or the arboreal gas model onto a hemi-supersphere sigma model: the correspondence with the 
sigma model holds exactly at the perturbative level, but breaks non-perturbatively. This is not 
expected to affect the physics in the broken symmetry phase where the field is slowly varying, 
and does not jump between the two solutions. The physics in this phase follows from the RG 
equation at one loop: 

^ = (5(g) = (n- 2)g 2 (12) 

with n = —1 = 1 — 2 for OSP(l/2) and g positive (i.e., a negative). We see that g flows to 
in the IR for small enough initial values, the theory remains massless, and the symmetry is 
spontaneously broken. 1 This is in complete agreement with the phase diagram of the Potts model 
in the limit Q~0, e K — 1 ~ ^ (see Fig. [5J). For the sigma model, the flow pattern holds until g 
hits a critical value, beyond which the model is massive in the IR again. On the square lattice 
on the other hand, the massless phase of the Potts model ends up with the antiferromagnetic 
transition point, and thus it is tempting to speculate that the universality class of the critical 
sigma model coincides with the critical antiferromagnetic Potts model. But this might be affected 
by the difference between the hemi-supersphere and the full one (the effect of global topological 
properties of manifolds on the critical points of sigma models does not seem to be completely 
understood in general. See for a recent review.) 

While the situation is not so familiar in two dimensions where most effort has been devoted to 
models based on unitary compact groups where the Mermin- Wagner- Coleman theorem prevents 
spontaneous symmetry breaking, it is well known in higher dimensions. There, one typically 

1 When g is negative (i.e., a positive), g flows to large values in the IR, where the model is massive, and the 
symmetry restored by the fluctuations. 
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e K > l-Q/4 : 
Massive region 

^ e K = l-Q/4 : 
r ~ - - . _ Critical point 

e K < l-Q/4 : " 
Massless region; 



flows to c=-2 



Figure 2: Phase diagram of the supersphere sigma model from the Q state Potts model on the 
square lattice. 



tackles exponents at the critical point in the O(N) model with ($ 2 ) 2 coupling by studying 
the broken symmetry phase, i.e., the non-linear sigma model, in d = 2 + e, calculating the (3 
function in this phase, and identifying the critical point with the first zero of the (3 function. The 
perturbative series thus obtained are only sensitive to the local structure of the group manifold, 
while the restoration of symmetry involves global properties. 

To discuss things further, we consider the geometrical interpretation of the full sigma model. 
It is slightly more convenient to start from the non-scaled fermions, and consider the equivalent 
of Eq. © with the constraint b(<fr) > removed. To each site i we assign an Ising variable ej = ±1 
which distinguishes between the two solutions (i.e., b(<p) > and b(<j)) < 0). One finds easily 



•^sigma 



Y[d(p(i)dr]i(i)dr} 2 (i) 5[cj) 2 (i) + 2ar]i(i)r] 2 (i) - 1] exp 



1 



-E(*(«"M(3)-i) 



(ij) 



Y[dr]i(i)dr] 2 (i) ^ exp 

{ei} 



<ij> 



(13) 



x exp 



(ij) 



where we have defined ej - 
and thus the existence of domain walls with e^e^ 



Hj^i m ij£j- 



For the region of broken symmetry, a is negative, 
— 1 extremely unfavorable for small \a\. If 
we think of the ej's as Ising variables (essentially the c- number parts of the field </>), we have a 
ferromagnetic Ising model coupled to the arboreal gas. A full geometrical interpretation is easily 
obtained by going over the steps of 0. The 77 terms again give an arboreal gas, but now each 
tree gets a more complicated weight 



a + 



i in tree 



(eiii - ma) - a (e^e,- - 1) 

(ij) in tree 



(14) 



so now trees getting near or across domain walls have their weight modified. 

As the lattice coupling is modified, a first sight analysis would say that generically in such 
a model, one will observe two transitions: an Ising transition where only the dynamics of the 
variables becomes critical, and an arboreal gas transition, where only the dynamics of the trees 
becomes critical. These transitions will occur for different values of a. Hence the critical point 
of the arboreal gas, i.e., the hemi-sphere lattice sigma model, will coincide in general with the 
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critical point of the coupled system [Ising model + arboreal gas], i.e., the full lattice sigma model. 2 
Only an extraordinary coincidence could make the Ising and arboreal transitions merge, resulting 
in a higher order phase transition, and a situation where the critical points of the hemi-sphere 
and full sigma model are different, the latter being of higher order, and less stable. 

This whole analysis may have to be taken with a grain of salt however. In particular, the 
arboreal gas is a model whose weights are extremely nondocal, and the modification brought by 
coupling to the Ising variables might have strange consequences. Also, it is after all not clear 
how much of the conventional wisdom derived from models with positive Boltzmann weights 
will extend to models such as ours, where the physics is determined by balance of positive and 
negative contributions in addition to the usual energy/entropy balance. At this stage, numerical 
study seems the best way to clarify the situation. 



3 Bet he equations at the antiferromagnetic point 

The special value a = —4 identified above corresponds more precisely to the end point Q — » of 
the antiferromagnetic critical line, which was identified by Baxter JJj as 

(e^ 1 + l) (e^ 2 + l) = 4 - Q (15) 

(here K\,K2 are the usual horizontal and vertical couplings). While the Potts model on the 
square lattice is not solvable for arbitrary couplings, it does become solvable on the line p5|). 
The Bethe ansatz equations were written in the original paper by Baxter, and read 

fsmh(X-i^ N TT sinhi(A-A / -2i 7 ) 

"11— TTT^ w , ^ ( 16 ) 



\sinh(A + i7)/ -^f sinh ~(A — A' + 2i^f) 

where N is the width of the cylinder in the transfer matrix (measured in terms of the Potts 
spins), the A are the usual spectral parameters determining the Bethe wave function, and we 
have parameterized \[Q = 2 cos 7, 7 G [0, |] . As discussed in details in p] for the case Q generic, 
the physics of the isotropic model K% = K2 can be studied more conveniently by considering a 
strongly anisotropic limit and the associated quantum hamiltonian for a chain of 2iV spins 1/2. 
The energy (in the relativistic scale) is then given by 

E _ K-2 1 y 2sin2 7 

A 

We expect that the ground state is obtained by filling up the lines Im A = 5 and Im A = — \. 
If we restrict to solutions of the type A = ^ + a and A = — ^ + f3, a system of two equations 
for the real parts a and j3 can then we written. It can also be extended to string configurations 
centered on either line. 



2 This is if the Ising transition has not yet occurred when the arboreal degrees of freedom become critical. The 
converse might happen, with the Ising variables becoming critical first, but this seems quite unlikely, as it would 
lead to the universality class of the supersphere sigma model being a simple Ising theory. Notice that at small 
negative a the Ising model is in the ordered phase, and the weight per unit length of a domain wall is e 2//a <C 1. 
As \a\ increases, this weight increases until at the critical value a = —4 it reaches e _1//2 = 0.606. Without coupling 
to the arboreal gas, the Ising model meanwhile would order when the weight per unit length of a domain wall is 
Y~^/2 = 0-414. So the Ising transition must be significantly hampered by the trees. 
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We now focus on the limit 7 — ► \. We set 7 = \ + e, and scale the spectral parameters with 
e. The resulting system of equations reads (we use the same notations for the rescaled spectral 
parameters) 

a - i \ N _ yr a - f3 - 2i 
a + i) ~ 11 a + (3 + 2i 
f3-i\ N Tr /3-a-2i 



and the energy 



13 + iJ a 13- a + 2i 
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4 Critical properties at the antiferromagnetic point 

4.1 First aspects 

We now observe that setting a = (3 gives Bethe equations similar to the antiferromagnetic XXX 
chain, while there is an additional factor of 2 in the definition of the energy. The sector symmetric 
in a and (3 roots will lead to twice the central charge and critical exponents of the XXX chain, 
which are the ones of a free boson with coupling constant g = 1, or an SU(2) WZW model at 
level one (a similar feature was encountered in [3]. This is more than a coincidence, as will be 
discussed below). 

This is not to say that all the excitations are of this type. Doubling the spectrum of the XXZ 
chain of course creates many gaps in the conformal towers, which have to be filled with other 
types of excitations in our system. But at least, we can derive quickly some results from this 
subset of the spectrum, which is now under control. 

For the electric part of the spectrum, we then get the following result 

c = 2 x ^1 - — = 2 - 12e 2 (20) 

The charge e can be interpreted as giving a special weight w to non-contractible loops. Following 
the detailed mappings of the Potts model on contours and vertex models, one finds w = 2cos7re. 
Since non-contractible loops come in pairs, we have e defined modulo an integer. The Q — > 
state Potts model itself corresponds to e = eo = |, for which c = — 1, while the effective central 
charge is c e ff = 2. 

Magnetic gaps read similarly 

2 2 
am m 

A = 2 x = — 21 

4 2 y J 

Here m is the spin of the equivalent XXX chain. Since one solution of the equivalent XXX 
system corresponds to a pair of solutions of our system, m = y — [number of pairs of solutions] 
and m = where S z is the spin in our initial system. Therefore 



A = ^ (22) 

where S z is an integer since we have an even number of spins. 

One can see here that parity effects will likely occur. If N is even, since we must also have an 
even number of upturned spins to make up pairs of Bethe roots in the "folding" of the equations, 
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the result strictly speaking will only apply to S z even, including the ground state. Numerical 
study shows that in fact it holds for S z odd as well. If N is odd on the other hand, since we must 
also have an even number of upturned spins to make up pairs of Bethe roots in the "folding" 
of the equations, the result strictly speaking will only apply to S z odd. Numerically, we have 
observed indeed that for S z even, including the ground state S z = 0, results for N odd obey a 
different pattern. This will be discussed later, and we restrict for now to N even. 

The electric charges in the S z = sector of the equivalent XXX chain are, for the untwisted 

2 2 

case, integer, and give rise to the following gaps A = 2 x |^ = Combining all the excitations, 

we thus get exponents from a single free boson at g = |, viz., A em = | [e ± , where e, S z 

are integers. 3 This is intuitively satisfactory, as one expects the "arrows degrees of freedom" in 
the model to be described by a free boson in general. This free boson of course contributes only 
c = 1 to the total central charge, so to get c = 2 we need to invoke the presence of another theory 
with c = 1. 

To proceed, we consider the dynamics of the hole excitations over the ground state where the 
two lines Im A = =t^ are filled. We introduce the functions (notations are the same as in 

a + in/2 i d 

en(«) = — ^1 a n (a) = --ln[e n (a)] (23) 

and the Fourier transform 

f(x) = / dae !a 7(«) (24) 



such that a n (x) = e n \ x \/ 2 . We take the logarithm of the Bethe equations [Eq. (|18|) ] and introduce 
densities of particles (p) and holes (p h ) per unit length. The result in Fourier transforms reads 

Pi + Pi = a 2 ~ 04 * P2 

h + p\ = a-2 - a 4 * Pi (25) 

and is of the form p + p h = p + Kp. Physical equations (that is, equations where the right hand 
side depends on density of excitations) meanwhile read 

h + ft = 2^sh^ + 1^4N % + 2 sinh 2x^ 2 

^ + &= 2^ X + 2A2-J" + T^r\& (26) 
-1 



Introducing the matrix Z = (1 — K) one finds 



-4W 2 sinh 2x 
-1 1 



(27) 



2 sinh 2x 



Now in ordinary models such as those with sl{n) symmetry, it is well known how to obtain 
analytically the spectrum of exponents by careful analysis of finite size scaling effects. The 
exponents are expressed as quadratic forms based on the matrix Z and its inverse |19| I20j. We 

3 In general, we will denote by A the (rescaled) gaps with respect to the central charge c = 2, and reserve the 
notation h for the (rescaled) gaps with respect to the true ground state of the antiferromagnetic Potts model. 
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shall assume a bit naively that the same essentially is true here as well. All the elements of Z 
diverge in the limit x = 0, but if we keep x small, Z reads 

Z ={ i?± ) ( 28 ) 

\ 4e 2 + 4e / 

which allows us to calculate the inverse matrix as e — > 

We can now let e — > to express the finite size scaling corrections and thus the gaps for N even 
as [20] 

x = -drfZ^Sn + d l Zd (30) 
4 

where 5n is the two-column vector ^ ^ , Jnj the change of solutions of the i th type. Similarly, 

d is the two-column vector ( j 1 ) , d the number of particles backscattered from the left to the 

\d 2 J 

right of the Fermi sea. Therefore we find 

A + A = -{8m + 5n 2 f + x (Sni - 5n 2 ) 2 + j(di + d 2 ) 2 + oo x (di - d 2 ) 2 (31) 

The symmetric part of the spectrum 8n\ = 5n 2 = 5n and d\ = d 2 = d gives h + h = (5n) 2 + d 2 . 
We have checked numerically that this formula extends to the case 5n half odd integer, so the 
total XXX part of the spectrum reads h + h = e 2 + e,m integers, or h = |(m + 2e) 2 
and h = |(m — 2e) 2 . This corresponds to a free boson at the coupling g = 1/2, i.e., symplectic 
fermions. It is useful to introduce now the determinants of the Laplacian on a torus with different 
types of boundary conditions, Det A a f,; the subindices denote the boundary conditions (P = 
periodic, A = antiperiodic, or F = free), with b referring to the space direction (of length N) 
and a to the imaginary time direction (of length M, with modular parameter q = e~ 2nM / N ) |21| : 

Det Aap = — T g(2n+l) 2 /8g(2m+l) 2 /8 

Det Aaa = -V^/s^/s 

' ' n,m 

Det Apa = — y(-l)V 2n)2/8 (-l) m g (2m)2/8 (32) 

where r\ = ? 1 ^ 24 ni°(l — Q n ) IS Dedekind's eta function. The generating function of XXX levels 
in the case iV even thus reproduces Det Aap + Det Apa + Det Aaa- A 

Combined with this symplectic fermion part, we have an infinite degeneracy associated with 
the modes d\ = d 2 . Presumably, this leads to a continuous spectrum of exponents on top of the 
symplectic fermion ones. Since on the other hand we know that the theory has c = — 1 (c c g = 2) 
while the symplectic fermion part identified so far has c = —2 (c e g = 1) , it is very tempting 
at this stage to propose that the missing part of the spectrum consists simply in a non-compact 
boson. We now give a more direct argument for this. 

4 Notice that since Det A PP = 0, we can write as well Det A A p = -4 V' q ^n+if /8^2m+i) 2 /s where the 
sum is over m and n of the same parity. 
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4.2 Lessons from free boundary conditions and modular transformations 

The partition function of the Q = Potts model on an annulus with free boundary conditions 
can be obtained in two steps, both based on previous works. The first step consists is expressing 
it as an alternate sum over generating functions of levels for type II representations of Utsl(2) 
(where \/Q = i + t _1 , and we use the notation t for the quantum deformation parameter to avoid 
confusion with q) 22 in the associated six-vertex or XXZ model (not to be confused with the 
XXX model which appears through its Bethe equations earlier in this paper as a subsector of 
the Potts model. The latter has t = 1 while the former has t = i). The second step consists in 
identifying the generating function Kj of levels for the type II representations. This step is not 
rigorous, and based on our understanding of the model at special points 7 = ~, n > 2, together 
with arguments of continuity in 7 |lUj . To make a long story short, one ends up with the following 
results. 

The expression of Kj as Q — » is 



o i0'+l)/2+l/24 
K — - 

" P 2( g ) 



00 00 

2n 2 +n(2j+l) _ r> „2(n+l/2) 2 +(n+l/2)(2j+l) 



1 + 2 q 2n +ni2j+1) ~2j2q 2 
n=l n=0 



(33) 



where P(q) = IIi^l ~~ 9 n )j an d j is an integer (spin) label, j = 0,1,2.... The astute reader 
might wonder why this expression does not satisfy for instance Kq — K\ + K2 — K3 . . . = 
which would naively result from the quantum group symmetry that mixes all spin and spin 
1 representations [22]. The reason is that the eigenvalues of the transfer matrix all go to zero 
as Q — > 0, and this transfer matrix cannot be rescaled in such a way that its elements remain 
all finite as Q — > 0. Alternatively, in a hamiltonian formulation, the ground state energy of the 
relativistic hamiltonian goes to —00 in the limit Q — > 0. In that limit, the Kj's are obtained after 
an additional renormalization, and do not have to satisfy the quantum group induced coincidences 
any longer. 

A little algebra shows that the generating function of levels (partition function in the confor- 
mal language) of the XXZ or six-vertex model on the annulus is 

Zfp- = £(2i + l)K 3 =I^j; q N ^' 2 (34) 
j=0 W N=0 

In the usual case of the critical Q = state Potts model, the equivalent of Kj would be 

o j0-l)/2+2/24 

*'= C 1 -*"") (35) 

(in that case one checks that Kq — K\ + K2 — if 3 . . . = indeed). The partition function of the 
vertex model would then be in that case 

2/24 00 

^vertex = 4 9 £ q N(N + l)/2 (36) 

It corresponds to the partition function of symplectic fermions on an annulus with antiperiodic 
boundary conditions in the time direction, det (—Aaf), and is as well a character of the c = 
-2 rational logarithmic theory [23j. The object Zf r F teyi written as Zfp ex = Z£f, tex X ± = 
det (—Aaf) x - can thus be interpreted as arising from a theory made of symplectic fermions 
and a non-compact boson, with partition function -. 
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The presence of the non-compact boson can be ascertained through a modular transformation. 
Writing Z P e p tex = 2 , we see that modular transformations in the usual Q = critical Potts 
model do close, since 



H-i/t) = yf^Fe 4 (-i/T) 

H-l/r) = \Z=*F0 2 (-1/t) 



t?(-1/t) = V^ V (t) (37) 
In the antiferromagnetic case meanwhile we now generate a continuous spectrum, according to 



vi-iM^J^^W (38) 

As for the arboreal gas or (the properly defined limit of) the Potts model itself, the partition 
function is 

^arboreal = £(-1)^. = « £ ( 2 AT + 1) (-1)^^+1)72 (3g) 

7V=0 

This can be written, using the Jacobi identity, as 

^arboreal = q l/Up(g) ( 40 ) 

The corresponding object in the usual critical Potts model (Zpp) would vanish identically. As 
before, Zpp oreal can be interpreted as the partition function of symplectic fermions (here, fermions 
periodic on the annulus, and the zero mode subtracted) times the contribution from a non- 
compact boson, -. 

We thus see that the study of the model with free boundary conditions greatly strengthens 
the claim that the continuum limit of the model is made up of a symplectic fermion theory and a 
non-compact boson. This of course agrees with c = — 1, c e s = 2. Let us discuss this identification 
further. 

4.3 The continuum limit of the Q = antiferromagnetic critical Potts model 

The action corresponding to this continuum limit reads (Boltzmann weight e~ s ) 

Jd 2 x [{d^f + 2d^ m d^ 2 ] (4i) 

with the boson (ft being non-compact. Now recall that the scaling limit of the arboreal gas, 
or the Q = Potts model in the Berker-Kadanoff phase, was identified with the OSP{\/2) 
supersphere sigma model, whose action reads exactly like Eq. (|41|) with the additional constraint 
that (p 2 + 2771772 = 1. To all appearances therefore, the critical point of the arboreal gas simply 
corresponds to a free OSP{\/2) field theory, where the constraint has been forgotten, and the 
symmetry is realized linearly. This interpretation requires some further discussion. The sigma 
model in the massless Berker-Kadanoff (spontaneously broken symmetry) phase has action 



S=-Jd 2 x [(d^f + 2^77^772] (42) 



g 

with g a negative quantity. It thus looks as if the natural action for the theory without constraint 
should be Eq. (|41j) with a negative sign. While this has no importance for the fermions (the 
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sign for them can be switched by a trivial relabeling) it is not so at all for the boson. While the 
model S = + / (dcp) 2 is well defined, the model with S = — J (d(j)) 2 is naively rather ill defined! It 
is possible to argue however, that since the theory is not compactified, this sign does not really 
matter. 

Start with the 'proper' free uncompactified boson, with action (we put numerical constants 
explicitly now) 

S=^Jd 2 xd^d^ (43) 

and thus propagator (4>(z, z)<fi(0)) = — 7— ln|;a;| 2 . The fact that <fi is not compactified allows all 
powers of <fi as fields, and hence by summation, all complex exponentials. However, since only 
normalizable states appear in the partition function (see |241 125j for further discussion on this 
point), only imaginary exponentials will show up there. This is compatible with the result of the 
functional integration 

Z = -t=-_=-_\ {qqf'da (44) 
V/mr r\r\ rjrj J-oo 

Hence in a transfer matrix study of such a model, we would see essentially a continuous spectrum 
above the central charge c = 1. Note how the presence of operators with arbitrarily negative 
dimension does not make the theory unstable, nor does it affect the ground state in the transfer 
matrix calculations. 

If we admit that the "full operator content" is built out of fields Q < y a+tb ) ( t > and derivatives, then 
a change of variables 4> ^ i9 should leave it invariant. This means that full operator content can 
be built as well out of e( a + lb ) 9 and derivatives, where now 9 has a propagator with the wrong 
sign, (8(z,z)6(0)) = — ^rln|z| 2 . Normalizable states now correspond to real exponentials, and 
their contribution to the partition function gives exactly the same result as Eq. (|44|) . which can 
therefore be considered as the partition function for the theory with action 

S = ~\\ d 2 xd^e (45) 

In other words, we define the functional integral in the case of Eq. (|45|) by analytic continuation, 
and since the set of allowed fields is invariant under this continuation, the partition function is 
the same in both cases. 

The fact that the critical point of the supersphere sigma model appears to be the free theory 
is a bit surprising. One would have expected instead the OSP{\/2) symmetry to be restored in a 
non-trivial way, and the model maybe to coincide with the critical 0(n = —1) model of Nienhuis 
|26| . discussed more thoroughly as an OSP(l/2) model in Ref. f3f- The latter model however has 
central charge c = — |, and definitely is very different from the one we obtain. This difference 
might be related to the difference between the hemi-supersphere and the full supersphere sigma 
model. It might also be that some of the simplifications in the definition of the solvable 0(n) 
critical model in (such as the definition of the Boltzmann weight) affect the physics in an 
unexpected way when n is negative, with, as a result, less universality than expected. 

5 Relation with OSP{2/2) theories 

5.1 osp(2/2) integrable spin chain. 

We now notice that the Bethe equations bear a strange similarity with equations written for a 
chain with si (2/1) symmetry. In particular, if we took a chain with alternating 3, 3 representations 
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(N of either type), the equations would read, in the appropriate grading [77115]. 



a- 2i\ N ^ra-(3-2i 



a + 2ij 1L a + f3 + 2i 
P-2i \ N = T r P-a-2i 

l3 + 2i) H/3-a + 2i 1 ' 

More generally, the Bethe equations for alternating representations build with n-supersymmetric 
tensors and the conjugate would look like Eq. (|46|) . but with the factor 2 replaced by 2n on 
the left-hand side (source terms). Our equations seem to correspond to taking "quarter spin" 
representations (in sl(2) inspired conventions, where the fundamental has spin 1/2). This in fact 
has a meaning, thanks to the work |28| where it was shown that the Bethe equations for a chain 
built with the four-dimensional typical representation (and L sites) read 

a - i(2b + 1)\ L _yra- (3 -2i 
a + i(2b + l)J ~~ ll a + + 2i 
a + j(2b-l)\ L _ u P-a-2i 

11 R— ,v-L 04 ^ > 



a-i(2b-l)J lL f3-a + 2i 

where b is the fermionic Dynkin number (denoted by 62 hi „28j. The fundamental of si (2/1) 
corresponds to b = —1, a degenerate case. The case 6 = meanwhile corresponds to the fun- 
damental representation of osp(2/2). Energy terms in this case can also be read from [2Hj and 
exactly match the term in the Q = Potts model case (after a rescaling of the Bethe parameters 
in |28j). We thus conclude that the Bethe equations and energy for the Q = antiferromagnetic 
Potts spin chain are identical with those for the integrable spin chain with osp(2/2) symmetry 
built on the four-dimensional fundamental representation [b = 0,j = 1/2]. 

Note that the length of the osp(2/2) system in these equations, denoted L in the Bethe 
equations, coincides with N, where the number of spin 1/2 representations in the spin chain is 
2N. Hence the size of the Hilbert spaces on which the two systems are acting are identical, and 
equal to 2 2N = 4 N . Since the Bethe equations also coincide, this is a very strong indication 
that the spectra of the two systems are identical, even though one has to be careful with such 
statements, as Bethe spectra are not always complete in the superalgebra case. 5 

5.2 Continuum limit of the osp(2/2) integrable spin chain. 

Now, as argued in |31| . the integrable model based on the four-dimensional fundamental repre- 
sentation of osp(2/2) sits inside the Goldstone phase of this model, and goes to the supersphere 
sigma model OSP(2/2)/OSP(l/2) in the scaling limit. The latter model can be described j3U] 
by the following action. Parameterize the S 1 ' 2 supersphere as 

fa = cos 4> (1 - 771772) 

(j) 2 = sin 4> (1 - 771772) (48) 
such that (f) 2 + 4>\ + 2771 772 = 1. The action is then 

S = - J d 2 x [(d^f (1 - 2771772) + 2^771^772 - 4771772^771^772] (49) 



5 Note that such a 'transmutation' from Bethe equations characteristic of a rank-one algebra to those of a 
rank- two algebra was met already in the sl(3) regime of the chain 12! )|. 
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where 4> is compactified, = + 2ir. The integrable chain corresponds to spontaneously broken 
symmetry, with the theory being free in the UV and massless in the IR, corresponding to a 
negative coupling constant g. Critical properties are described by the g —* limit of Eq. (|49|) . A 
rescaling and relabeling brings Eq. Q49JI into the form 

s = Jd 2 x [- {d^f (i + 2\ g \r hm ) + w^d^ + Mglmm^md^} (50) 

with now <b = 6 -\ — 7=. As o — > 0, we thus obtain a symplectic fermion and a non-compact 

boson, with the wrong sign of the metric. This is exactly the low-energy theory we had identified 
earlier for the Q = antiferromagnetic Potts model, thus confirming our picture. 

Note that the presence of a continuous spectrum in this osp(2/2) version of the problem has 
a physical origin in the fact that the symmetry is spontaneously broken, and thus correlation 
functions of order parameters have no algebraic decay (though have logarithmic behaviour). 
This feature was verified numerically in |31| where it was found that the scaling dimensions of 
p leg polymer operators extracted through transfer matrix calculations vanish inside the entire 
Goldstone phase indeed. 

As a final note we remark that the Bethe equations for the osp{2/2) chain have been most 
often |32M33| written in a different grading from the one we used. It is also interesting to compare 
the analysis of the low energy physics from the more standard point of view. 

The usual Bethe equations are written as 

a — i\ L -j-jCt — /3 — 2z 



with the energy 



a + i/ LL a — (3 + 2i 

yr (3-a-2i _ (3-f3'-Ai 

i - i -(3-a + 2i (3 - f3' + Ai { ' 

£ = -V-^ (52) 
7T ^ 1 + a 2 v ' 

(note the change of sign compared to the other Bethe ansatz). The relevant solutions for the low 
energy physics are: a two-strings centered on real /3's (denoted then by y), and real a's (which 
are not in complexes, denoted then by z). The corresponding Bethe equations are 

y-3iy + i \ L _ y. y-y'-Ai y. y-z-2i 
y + 3iy - ij y — y' + 4i y — z + 2i 

z — i\ y-r z — y — 2i 



n+s^ pa) 



z + ij L y L z-y + 2i 

The ground state is filled with the complexes. Going to the continuum limit and using Fourier 
transform gives the following physical equations, if we call p the density of complexes, a the 
density of real particles 

1 q-\ 2x \ „ h 1 

^ ^ 2 cosh x 2sinh2|x|^ 2sinh2|x| 

1 e" 2 ! 1 '! 1 

a + a h = , , - a + ^^p h (54) 

2 cosh x 2smh2 \x\ 2sinh2x 
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These equations show that the chain has a relativistic limit (a feature peculiar to the value b = 
of the fundamental representation Moreover, the density equations coincide with those for 

the Q — * antiferromagnetic Potts model discussed above: even though we have started from two 
different gradings, we end up with identical physical equations, a very strong indication that our 
approach is consistent [31] (note in particular that the quantum numbers of the excitations, being 
expressed naturally in terms of complexes and real holes, have different expressions in terms of 
the bare quantum numbers) and the results correct. 

To finish the analysis of the osp{2/2) system, the most useful quantity one can consider is 
the 'equivalent spin' defined by forming the combination m = b — 3S Z . The four states of the 
fundamental representation b = 0, S z = ±1/2; S z = 0,6 = ±1/2 correspond respectively to the 
values m = 3/2,1/2,-1/2,-3/2. On the other hand, according to the detailed lattice model 
derivation of the Bethe ansatz equations [33] EH] , 

m = b - 3S Z = #/3 - #a - - 

The spectrum is expressed in terms of sums and differences of 5n^ and 5n r , where the subscripts 
c, r indicate complexes and real particles respectively. One has 

5ric — 5n r = 5np — 5n a 

dn^ + 5n r = bnp — 35n a (55) 

We thus see that it is possible to change the value of the pseudospin m without changing 5n^ + 
5n r , i.e., by exciting only the flat part of the spectrum. In other words, in the geometrical 
interpretation, the gaps for p leg operators are all expected to vanish. This result agrees with 
predictions based on the Goldstone analysis of the OSP{2/2) loop model, as well as numerical 
calculations [5J, and we take the consistency of the different approaches as a justification of our 
use of the dressed charge formalism in these problems. 

5.3 Explicit mapping between the two models 

Even though we have found that the Bethe equations for the two models are identical, we do not 
exactly know how to map them onto each other. Clearly, one must combine pairs of neighbouring 
spin 1 /2 representations of the vertex model to make up four-dimensional representations of the 
osp{2/2) algebra, but we are not sure how to implement this. Part of the difficulty comes from 
the highly singular nature of the Potts transfer matrix as Q — > 0. 

It is interesting to notice however that as t — » i, the tensor product of two fundamental 
representations of Utsl(2), which decomposes as the sum of a one and a three dimensional repre- 
sentations like for t = 1 when t is generic, becomes indecomposable, and has the structure shown 
in Fig. 01 where |u a ) = e i7r / 4 | + -) + e" i7r / 4 | - +) and \u 2 ) = e~ in / 4 \ + -) + e i7r / 4 | - +}, and the 
arrows indicate the action of the Utsl{2) generators: 

S~ = t° z ® a' + a' (8) t~ aZ 

S + = t aZ ®o + + a + ® t-° z (56) 
It turns out that, if one introduces the other generators 

£T = r° z <g> cr" + a~ ® t aZ 

S + = t~ aZ ® <7+ + <j+ ® t° z (57) 
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Figure 3: Indecomposable representation 1/2® 2 for t = i. 





Figure 4: Mapping between forests and self-intersecting loops. Thin solid (resp. dashed) lines 
represent the lattice C (resp. the dual lattice £ D ). All edges of CP penetrating the boundary of C 
are understood to join at a unique dual site situated on the infinite, exterior face of C. Thick solid 
lines depict the forest (here with four component trees); thick dashed lines depict the (unique) 
dual tree. The grey (red online) solid line is a self-intersecting loop on the surrounding lattice 
(see text). It bounces off occupied edges (thick linestyle) on C and on C D , and self- intersects 
whenever a pair of intersecting direct and dual edges are both non-occupied. 



at t = i, together with the limits lim^j \ +t -i ? one obtains the four-dimensional representation 
of osp(2/2). This suggests a mechanism by which the symmetry of the lattice model must be 
enhanced at the critical point. One can in fact show by brute force that by concatenating pairs 
of spins 1 /2 into four degrees of freedom and carefully taking the Q — > limit of the Boltzmann 
weights of the heterogeneous 6-vertex model equivalent to the critical antiferromagnetic Q state 
Potts model, the R matrix of the osp{2/2) integrable vertex model in the fundamental ([b = 
0,j = 1/2]) representation is indeed obtained. This will be discussed elsewhere. 

Note of course that the osp(2/2) integrable model can be reformulated as a loop model. In fact 
the whole Goldstone phase has been argued to be realized physically by taking a model made of 
a single dense loop where crossings are allowed. This can be obtained (and studied numerically) 
for instance by taking the usual [HE! covering of the surrounding lattice of the square lattice with 
a single self avoiding loop, and then allowing self-crossings. 

It is well known that in the case where self-crossings are not allowed, the loop the hull of a 
spanning tree, hence giving rise to the well-known close relationship between the critical theory 
of dense polymers and that of spanning trees, both being related to symplectic fermions |12| I13j. 
This relationship can to some extent be generalized to the arboreal gas and the loop with self- 
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crossings (see Fig. |1J). Indeed, draw on the lattice C a spanning forest. Draw on the dual lattice 
C D (which in particular contains a single site outside C) a spanning tree. The resulting contour 
on the surrounding lattice is a single loop with a number of self-intersections equal to the number 
of trees on C plus one! So we have a sort of correspondence between a self-crossing loop and an 
arboreal gas, i.e., between a system with osp(2/2) symmetry and the Q —* antiferromagnetic 
Potts model. However, this is not a very satisfactory correspondence. In particular, because 
the tree on the dual lattice plays a crucial role in determining the loop configuration, there is 
no one to one correspondence between loop and forest. Also, giving a negative fugacity to trees 
corresponds to a negative weight for self-crossings of the loop, while the existence of the Goldstone 
phase was argued in jHJ for positive weights only. 

6 Conclusion 

The most interesting conclusion is probably the emergence of a non-compact boson in the con- 
tinuum limit of the critical antiferromagnetic Potts model at Q = 0. Although it is difficult 
to demonstrate the existence of such a non-compact object based on any single analysis, the 
convergence of Bethe ansatz calculations, numerical calculations, and symmetry considerations, 
leaves very little doubt that the conclusion is correct. This is, to our knowledge (together with 
our earlier work on Goldstone phases which as we saw is intimately related), the first time 
that a non-compact object appears in the continuum limit of a lattice model with a finite number 
of degrees of freedom (but 'unphysicaF Boltzmann weights). In a subsequent work, we will see 
that this feature extends to the antiferromagnetic Potts model away from Q = j^j. 

We have also seen how this non-compact boson is part of a free OSP(l/2) theory, which 
appears here as a critical point of the super (hemi-)sphere sigma model. The main question that 
remains unanswered in this study of the physics of models with supergroup symmetry is what 
happens to the 'other' (Nienhuis) critical 0(n = —1) model |26| . and how this is related with 
a possible difference between the critical points in the hemi and full supersphere sigma models. 
We hope to get back to these questions soon. 

To finish, we can clarify several remaining points. The first is the issue of what happens 
when the size of the lattice system N is odd. As will be discussed in details for the general 
antiferromagnetic Potts model in the spectrum in that case does not show any signs of 
non-compactness because the bosonic field </> appears to be twisted (and thus its properties 
independent of any compactification radius). By OSP{\/2) symmetry, the fermionic field is also 
twisted, so the effective central charge in this sector is 

c efr = c - 24/4 - 24C = -l-| + 3 = | (58) 

in agreement with numerical calculations [HUH]- 

The second point is that numerical [H] as well as analytical arguments |5] suggest that the 
critical point is, in the Potts model point of view, approached from the Berker-Kadanoff phase 
(the broken symmetry phase) with a singularity in the free energy of the form 

/ oc (a* — a) ln(a* — a) as a f a*. (59) 

(We recall that a* = —4 for the square lattice.) Using hyperscaling, this corresponds to a 
perturbation of the critical point by a field of vanishing dimension. There are of course plenty of 
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a 


L = 6 


L = 8 


L = 10 


L = 12 


L = 14 


xy. 





0.0000 


0.0000 


0.0000 


0.0000 


0.0000 




-1 


-0.0859 


-0.0792 


-0.0747 


-0.0714 






-2 


-0.1092 


-0.0992 


-0.0925 


-0.0877 






-3 


-0.1168 


-0.1079 


-0.1003 


-0.0947 






-4 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 







0.7492 


0.7498 


0.7499 


0.7500 






-1 


0.5823 


0.5975 


0.6075 


0.6146 






-2 


0.5216 


0.5492 


0.5663 


0.5781 






-3 


0.4936 


0.5253 


0.5464 


0.5609 






-4 


0.5466 


0.5691 


0.5830 


0.5927 




x 3 : 





1.9671 


1.9894 


1.9958 


1.9980 


1.9989 




-1 




1.6959 


1.7213 


1.7384 






-2 


1.4930 


1.5727 


1.6201 


1.6514 






-3 


1.4050 


1.5062 


1.5667 


1.6065 






-4 


1.5268 


1.6157 


1.6621 


1.6907 


1.7105 



Table 1: Scaling exponents x^ = 2h^ describing the probability P(r) ~ r~ 2xk that k distinct 
trees come close in points and r. The data is obtained by fitting transfer matrix results in 
the random-cluster picture, on periodic square-lattice cylinders of widths L and L — 2. Note the 
strong finite-size dependence for tree fugacity a^0. For a = the exact result is x^ = (k 2 — 1)/4 



such fields in the non-compact boson theory. We speculate that the scaling limit near and below 
a* could be described by the simple Landau-Ginzburg theory 

S = J d 2 x [{d^f + 2d^ m d^ 2 ] + (a* - a) (f + 2 mm - l) (60) 

for which the logarithm in / [see Eq. (|59)) ] follows as well. 

The last point concerns the physical interpretation of our results. This of course has to be 
taken with a bit of care, as our model involves negative Boltzmann weights. But the image we 
have is that in the region a* < a < 0, the partition function is dominated by configurations where 
a single infinite tree covers a finite fraction of the lattice — this fraction is unity as a — > 0, and 
remains finite in the whole phase a* < a < 0, the rest of the lattice being covered by trees of 
finite size. Tree exponents are the same as in the limit a — ► 0, as can be checked numerically 
in Table 1. In particular, the fact that the one-tree exponent remains equal to zero guarantees 
that there is an infinite tree, and that its fractal dimension is equal to two. At a = a* , we have 
coexistence of two phases, one with an infinite tree of finite density, the other with only finite 
trees. As a = a* is crossed, the derivative of the free energy with respect to a is discontinuous jSj, 
implying that the average number of trees per lattice site, and thus the average size of finite trees, 
experiences a discontinuity. In the phase a < a* (like in the phase a > 0), there are only finite 
trees, of typical linear size given by the finite correlation length. The point a = a* is actually a 
first-order critical point, that is a point where a first and second order phase transition coincide 
(Hj. At such a point, some correlations still decay algebraically, while some of the thermodynamic 
properties have discontinuities. The latter translate into operators of vanishing dimensions, for 
which we expect an infinite multiplicity. The simplest formalism to encompass this multiplicity 
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is a non-compact boson, in agreement with what was found previously. 

We conclude with a side remark. The fermionic representation of the arboreal gas model, 
instead of being used to argue for the presence of an OSP(l/2) symmetry, can be used to get yet 
another representation, this time in terms of a loop model. For this, go back to the initial path 
integral © written in terms of fermions rji , 772 only 



exp 



(^^E^W^W - ^2(vi(i)V2ti) +m(j)V2(i)) — * »7i (0*72 (*)»7i 0")»72 0/ ) 



(61) 



Simply expand the exponent, and contract the fermions. This is like the standard high tempera- 
ture expansion for the Nienhuis 0(n) model . It gives rise to a sum over non- intersecting loops, 
with weight —2 per loop (like in the 0(n = —2) model), plus non-overlapping monomers with 
weight 4 + a, and dimers (which can be considered as degenerate length-two loops) with weight 
—a. The critical point a = —4 means monomers disappear, dimers have weight 4. We are not 
sure why this model should be critical — note the dimers seem to encourage swelling of the loops, 
instead of self attraction as one would expect to get to a tricritical point — and this direction 
certainly deserves further study. Note that a critical point with c = — 1 has been identified in a 
modified version of the 0(n = —2) model in Ref. |36| . 

Acknowledgments: This work was supported by the Department of Energy (HS). We thank 
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A Bet he equations 

The standard integrable model based on the s/(2/l) = osp{2/2) superalgebra is the Perk Schultz 
model. It uses the fundamental representation of dimension 3, and leads to the Bethe ansatz 
equations (for a review see jSZ]) 

'a -2i\ N j*a-0-2i 



a + 2ij 11 a- (3 + 2i 
-pj- 0-a-2i = -pj- - 0' - 4j 

11 0-a + 2i 11 0-/3' + ii [ ' 

Algebraically, these equations correspond to treating the fundamental representation of s/(2/l) as 
a non- typical fermionic representation with parameters [b,j] = [—1/2, 1/2], of dimension 3 instead 
of 4. One could as well consider a model based on higher spin representations [b,j] = 
obtained by taking the fully supersymmetric component in the product of n = 2j fundamental 
representations. The equations would read as Eq. ()62j) but on the left-hand side the factor 2 
would be replaced by 2n. 

More generally, for the four-dimensional typical representations [b,j = 1/2], the Bethe equa- 
tions generalizing Eq. (|62|) read [SHI OH ED] 

'a + i(2b-l)\ N rja- (3-2i 



a - i(2b -1)7 II a - + 2* 
0-a-2i _ p-£-U 
ll/3- a + 2z -^0-0> + Ai {M) 

In particular, the case b = is similar to Eq. ()62|) but on the left-hand side the factor 2 is replaced 
by 1 corresponding roughly to taking a 'quarter spin'. 
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The second type of Bethe ansatz for the Perk Schultz model is obtained when associating 
the fundamental representation to the bosonic root of the distinguished Dynkin diagram instead. 
One obtains then the so called Sutherland's equations 

_yr a-(3-2i 
i - i -a-(3 + 2i 

P-2i \ N T[ (3-a-2i _ n /?-/?'- 4z 

(3 + 2i) i - i -(3-a + 2i 11 + { ' 

Finally, it is also possible to write a Bethe ansatz based on the alternative Dynkin diagram, where 
both roots are fermionic 

1 = jt a- (3-2i 
H a- (3 + 2i 
(3-2i\ N 1 - f (3-a-2i 



f3 + 2l J n /3 _ a + 2i ( 65 ) 
In this ansatz, equations for a mix of the fundamental and its conjugate are particularly elegant 

a - 2i \ N _ -pj- a - (3 - 2i 
a + 2i) ~ H a - (3 + 2* 

JVM) =^h-a + 2 i (66) 

As for the model based on the four-dimensional representation, its equations are given in the text 
as Eq. Q47|) and correspond once again to taking a 'quarter spin'. 

References 

[1] M. Zirnbauer, Conformal field theory of the integer quantum Hall plateau transition, 
|hep-th/9 905054 

[2] I. A. Gruzberg, A. W. W. Ludwig and N. Read, Phys. Rev. Lett. 82, 4524 (1999). 
[3] N. Read and H. Saleur, Nucl. Phys. B 613, 409 (2001). 

[4] V. Gurarie and A. W. W. Ludwig, Conformal field theory at central charge c = and 
two-dimensional critical systems with quenched disorder, |hep-th/0409105] 

[5] F. Essler, H. Frahm and H. Saleur, The continuum limit of the integrable si (2/1) spin chain 
with alternating 3,3 representations, cond-mat/0501197 

[6] R. M. Gade, J. Phys. A 32, 7071 (1999). 

[7] S. Caracciolo, J. L. Jacobsen, H. Saleur, A. D. Sokal and A. Sportiello, Phys. Rev. Lett. 93, 
080601 (2004). 

[8] J. L. Jacobsen, J. Salas and A. D. Sokal, Spanning forests and the q-state Potts model in the 
limit q -> 0, |cond-m at/0401026 . 



21 



[9] J. L. Jacobsen and H. Saleur, in preparation. 
[10] H. Saleur, Nucl. Phys. B 360, 219 (1991). 
[11] R. J. Baxter, J. Phys. C 6, L445 (1973). 
[12] B. Duplantier and H. Saleur, Nucl. Phys. B 290, 291 (1987). 
[13] B. Duplantier and F. David, J. Stat. Phys. 51, 327 (1988). 
[14] H. Saleur, Nucl. Phys. B 382, 486 (1992). 

[15] R. J. Baxter, H. N. V. Temperley and S. E. Ashley, Proc. Roy. Soc. London A 358, 535 
(1978). 

[16] D. Mouhanna, Autour des systemes classiques magnetiques frustres, these d'habilitation, 
Universite Pierre et Marie Curie, Paris VI. 

[17] R. J. Baxter, Proc. Roy. Soc. London A 383, 43 (1982). 

[18] H. Saleur, Nucl. Phys. B 578, 552 (2000). 

[19] J. Suzuki, J. Phys. A 21, L1175 (1988). 

[20] H. de Vega and E. Lopes, Nucl. Phys. B 362, 261 (1991). 

[21] P. Di Francesco, P. Mathieu and D. Senechal, "Conformal field theory" (Springer Verlag, 
New York, 1997). 

[22] V. Pasquier and H. Saleur, Nucl. Phys. B 330, 523 (1990). 

[23] A. Bredthauer and M. Flohr, Nucl. Phys. B 639, 450 (2002). 

[24] K. Gawedzki, Non-compact WZW conformal field theories, hep-th/911076. 

[25] H. Saleur, in preparation. 

[26] B. Nienhuis, J. Stat. Phys. 34, 731 (1984). 

[27] J. Links and A. Foerster, J. Phys. A 32, 147 (1999). 

[28] M. Pfannmiiller and H. Frahm, J. Phys. A 30, L543 (1997). 

[29] L. Mezincescu, R. Nepomechie, P. K. Townsend and A. Tsvelik, Nucl. Phys. B 406, 681 
(1993). 

[30] B. Wehefritz-Kaufmann and H. Saleur, Nucl. Phys. B 628, 407 (2002). 

[31] J. L. Jacobsen, N. Read and H. Saleur, Phys. Rev. Lett. 90, 090601 (2003). 

[32] M. J. Martins and P. B. Ramos, J. Phys. A 28, L525 (1995). 

[33] P. B. Ramos and M. J. Martins, Nucl. Phys. B 500, 579 (1997). 

[34] F. H. L. Essler, V. E. Korepin and K. Schoutens, Int. J. Mod. Phys. B 8, 3205 (1994). 



22 



[35] R. J. Baxter, Exactly solved models in statistical mechanics (Academic Press, London, 1982). 

[36] H. W. Blote and B. Nienhuis, J. Phys. A 22, 1415 (1989). 

[37] F. Essler and V. Korepin, Phys. Rev. B 46, 9147 (1992). 

[38] Z. Maassarani, J. Phys. A 28, 1305 (1995). 

[39] M. P. Pfannmiiller and H. Frahm, Nucl. Phys. B 479, 575 (1996). 

[40] A. J. Bracken, M. D. Gould, J. R. Links and Y. Z. Zhang, Phys. Rev. Lett. 74, 2768 (1995). 



23 



